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The generalized purity is employed for investigating the process of coherence loss and derealiza- 
tion of the Q-function in the Bloch sphere of a two-mode Bose-Einstein condensate in a symmetrical 
double well with cross-collision. Quantum phase transition of the model is signaled by the gen- 
eralized purity as a function of an appropriate parameter of the Hamiltonian and the number of 
particles (N). A power law dependence of the critical parameter with N is derived. 
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Recently it has been pointed out that any bipartite 
and multipartite entanglement measure can signalize the 
presence of a quantum phase transition (QPT) [TJ [2j [3] 
in many particle systems. Related to that, a subsystem- 
independent generalization of entanglement has been in- 
troduced based on coherent states and convex sets char- 
acterizing the unentangled pure states as coherent states 
of a chosen Lie algebra [4]. Such a notion of entangle- 
ment denned relative to a distinguished subspace of ob- 
servables is pointed as particularly useful for classifying 
multipartite entanglement and thus for QPT character- 
ization. The generalized purity (GP) of the state rela- 
tive to a certain distinguished subset of observables form- 
ing a local Lie algebra is directly related to the Meyer- 
Wallach measure [5 j, and whenever a specific subsystem 
can be associated to the subset of observables the usual 
entanglement notion is recovered. On the other hand 
it has been demonstrated that the model describing a 
two mode Bose-Einstein condensate (BEC) in a sym- 
metric double well (BECSDW) [6 j, with cross-collisional 
terms [7J [8] presents interesting dynamical regimes: The 
macroscopic self-trapping (MST) and Josephson oscilla- 
tions (JO) of population [6j [8j [9], both experimentally 
observed [10 . Recently a discussion of the transition 
from one regime to the other has been presented from 
the point of view of critical phenomena, as a continuous 
quantum phase transition problem in terms of the usual 
subsystem entropy [TTJ [12J [13] . 

In this paper we apply the concept of GP in a 
BECSDW. Our purpose is twofold: (i) We employ the 
GP relative to a chosen set of observables to quantify the 
quality of the semiclassical approach used in most of the 
treatments of BECSDW; (ii) We characterize the quan- 
tum critical phenomena occurring in this model with the 
subsystem independent measure of quantum correlations 
[4]. The two mode approximated BECSDW has been well 
studied in the literature (see e.g. [6j [9] ) as a model 
presenting a nonlinear self-trapping phenomena. More 
recently, the importance of the cross-collisional terms for 
large number of particles in the condensate (TV ^> 1) has 
been noticed [7] and explored semiclassically with a time 
dependent variational principle (TDVP) based on coher- 



ent states |S]. For a fixed number of condensed particles 
TV, in order to explore the natural group structure of the 
model, we conveniently adopt the Schwinger's pseudo- 
spin operators defined in terms of the creation and an- 
nihilation boson operators d\_ , d± on the approximated 
localized states \u±) [6j [8]: J x = — d\d+)/2, 

J y = z(dld + - d\_d-)/2, and J z = {d\ d_ + d ] _d+)/2, 
where J = N/2. In that form the two- mode BEC Hamil- 
tonian writes as 



H = 2 



2A(N - 1) 
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where ft is the tunneling parameter, k is the self-collision 
parameter of the condensate which is much larger than 
the so called cross-collision terms n = fte 2 , A = ^e2, 
with e = (u+\u-). ft' ee 2[2A(7V - 1) + ft/2] is an effec- 
tive tunneling parameter dependent on N. The natural 
associated algebra of the model is su(2). In this form the 
Hamiltonian dTl is a realization of the Lipkin-Meshkov- 
Glick model [14], whose ground state entanglement has 
been investigated recently [13] in the 77 = limit . 

By means of a semiclassical method exploring SU(2) 
coherent states we have shown that for a sufficiently 
large N even a small amount of cross-collisional rate can 
change the dynamical regime of oscillations of the con- 
densate from MST to JO [8]. The quantum and semi- 
classical MST and JO dynamics of the BECSDW are 
known to be qualitatively very similar [6j [8] for large 
number of particles (N ^> 1), except for the presence of 
collapses and revivals in the quantum time evolutions for 
the mean values of (J x )(t). This breaking of quantum- 
classical correspondence is due to decoherence [18] intro- 
duced by quantum fluctuations, which drives the state of 
the system away from a coherent state. When we treat 
the semiclassical dynamics employing a TDVP [T5] . we 
restrict the evolution of the state to a nonlinear subspace 
constituted only by the coherent states of SU (2) [8] . Such 
an evolution is exact for an initially coherent state only in 
the macroscopic limit of TV — » 00 [T9] , Also, in the TDVP 
approach we force the system wave function to be always 
localized in the phase space. From this point of view, 
the coherent states are the closest to the classical ones, 
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which are points in the classical phase space, and con- 
verges to them in the limit N — > oo. The derealization 
of the wave function (DWF) which goes along with the 
decoherence, is the responsible for the quantitative dis- 
agreement between the two dynamics [8] , when collapses 
and revivals of expectation values of relevant observables 
happens with the DWF and consequent self-interferences 
occasioning the appearance of superposition states. 

For a more complete analysis of the decoherence or of 
the quality of the semiclassical approximation, it becomes 
fundamental to have a good quantitative measure of the 
'distance' of a given state to the subspace of the coherent 
states that give birth to the classical phase space. When 
the complete dynamics of the system is restricted to a 
space that carries an irreducible representation of SU (2) 
(here this condition is satisfied due to the particle number 
conservation) , there is a simple measure called GP of the 
algebra V su {2) [UHS], defined for su(2) as 

v su{2) m = ji E ww?- ( 2 ) 

k=x,y,z 

It is a good measure of decoherence, among other rea- 
sons, because it is invariant under a transformation 
of the group SU(2) on the state V su {2){\iji)) = 

V su ( 2 )(U\il))), V U e SU(2). Thus, all the states con- 
nected by a SU(2) transformation possesses the same 
purity. However, the most interesting property of V su {2) 
- concerning the purpose of quantitatively compare the 
correspondence between the semiclassical and the exact 
quantum dynamics - is the fact that this measure has 
its maximum value at one, if and only if, the state is 
the coherent state closest to the classical state: | #,(/>) = 
R{p,<t>)\J,-J) = e-^( J - sin ^- J y cos ^)|J,-J). As soon as 
such a state moves away from the coherence, becoming 
delocalized in the phase space, the GP decreases mono- 
tonically to zero. Remark that the GP only has such 
properties clearly defined for pure states and, only in this 
case, it is a measure of existing quantum correlations of 
the state on the classical phase space. 

We choose \J,J) X = \0 = f,0 = 0) where N 100 
particles are in the same well, Q = 1 and k = For 
this set of parameters it is known that choosing n = ^ 
the state is in the self-trapping region of the phase space, 
whereas for n = it is outside the self-trapping region 
and thus in the JO regime [8]. In Fig.Q we plot the GP, 
where the MST regime is in dashed line, and solid one 
for JO regime. In the MST regime the GP quickly drops 
down from 1 stabilizing at 0.9 for Qt ~ 10, indicating that 
the dynamics takes the state away from the subspace of 
CS. This plateau of purity coincides with the collapse 
region of population dynamics [8 . Note that close to 
Qt = 30 the purity shows small oscillations and in the 
region close to Qt = 60 the purity increases again until 
it practically recovers the value 1. At this instant the re- 
coherence of the state happens, being responsible for the 




°0 10 20 30 40 50 60 70 SO 90 100 

at 



FIG. 1: (color online) Time evolution of the GP V su (2) for 
MSP and JO dynamical regimes of the two mode BEC model. 
The dashed curve represents the MST regime with k = ^ 
and n — ^ , whereas the solid curve represents the regime of 
JO with k = ^ e n = For both cases the initial state of 
the system is the CS | J, J) x = \6 = § , = 0), where N 100 
particles are initially in the same well. 



revival of the oscillations of the population dynamics [8] . 
The agreement between the quantum and semiclassical 
results is due to the oscillation of the mean value of the 
generator J x in the MST, since the purity depends on the 
normalized square of such mean value. The oscillations 
of (J x ) around a non-zero value close to J keeps also the 
value of V su {2) close to its maximum possible value. In 
the JO regime for n = ^ and rj = the purity decays 
rapidly with the time, but in this case the decoherence is 
much stronger, and the purity reaches much lower values, 
not recovering to 1. Hence, this regime presents lower 
quantitative agreement between the quantum and clas- 
sical evolutions, when compared to the MST regime. In 
this dynamical regime, the system does not recover high 
values of coherence. At a time close to Qt = 250 the state 
reaches its maximum re-coherence, but the GP value 
is still lower than 0.4. This result is expected because 
the classical orbit delocalizes much more on the Bloch 
sphere for this regime, and this entails a correspondingly 
large derealization of the semiclassical Q-function on the 
sphere, Q(0, 0) = (0, (j)\p\0, (/>), with p = \^){^\, during 
its evolution [HI [17]. The larger the region traveled by 
the trajectory in the phase space, larger the broaden- 
ing of the distribution and smaller the coherence left on 
the state. Therefore, when the Hamiltonian has nonlin- 
ear terms in the generators of the dynamical group, the 
semi-classical approximation has better quantitative ac- 
cordance with the exact quantum results (for finite N) 
when the classical orbits sweep smaller 'volume' in the 
phase space. Since BECSDW model is integrable, we 
cannot analyze the decoherence due to chaotic trajecto- 
ries. However, our results indicate that the semiclassical 
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method has lower validity for this type of trajectory that 
is less localized in the phase space. 

Now we can take advantage of such qualities of the GP 
to characterize the quantum phase transition (QPT) [3]. 
The QPT is connected to a non-analyticity of the energy 
of the fundamental state of the system, when it is taken 
as a function of some real continuum parameter of the 
Hamiltonian [1 at zero temperature, in the thermody- 
namic limit N — > oo. Generally speaking, the energy of 
the ground state in a finite system is an analytic func- 
tion of any parameter of the Hamiltonian and only shows 
non-analyticity when N —> oo, which then corresponds 
to an avoided level crossing. However, even when we are 
not allowed to take effectively such a limit, we still can 
observe the scaling of the properties of the system for 
increasing TV and infer about the occurrence of the QPT 
in the thermodynamic limit. 

The BECSDW suffers a sudden change in its dynam- 
ics when k c = in the limit of no cross-collision terms 
and N ^> 1. With ^, ft > 0, such a transition of regime 
(JO to MST) does not occur at the ground state, but at 
the largest energy state for the value of parameter corre- 
sponding to the bifurcation in the classical phase space, 
which causes the appearance of a separatrix of motion. 
Strictly speaking, the QPT is only characterized in the 
limit N —> oo, thus our transition of dynamical regime 
(even if it occurred at the ground state) would only be 
considered as a bona fide QPT in the exact classical limit. 
In Fig.Q, we have the Q-functions for the eigenstate of 
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FIG. 2: (Color online) Q-function for the largest energy state 
of the H spectrum for several value of parameters k and n. 
Considering the number of particles N = 100 and Q = 1, we 
have the following values for the collision rates: (a) k = n = 0; 
(b) K = ^,77 = 0; (c) K = §, V = and (d) K =^, V = ^. 

largest energy of the spectra of H for several values of pa- 



rameters ft and rj. Fig.Qa just shows the coherent state 
\6 = 7r, <p) = | J, J) z , which is the maximum energy state 
of the non-interacting case ft = n = 0. In the absence 
of collisions, the most energetic eigenstate corresponds 
to the most localized state in the phase space, such that 
^sn(2) = 1- Increasing ft, but still not considering the 
cross-collision terms, the Q-function broadens along the 
x-axis, and consequently we expect the decreasing of the 
GP. At ft — 

shown in Fig. ^ 6, the state is greatly 
broadened, but still does not show a bifurcation; namely, 
the formation of two maxima in its distribution. This 
behavior is expected, since for finite N, the quantum 
transition parameter k%(N) is slightly different from the 
value of transition ft c of the classical limit. But for in- 
creased ft, such as in Fig.([2])c, for ft = ^ we see two 
maxima far apart along the x-axis as a signature of the 
bifurcation. The increase of the cross-collision has the 
opposite effect. For ft = r\ = the two peaks of the 
Q-function become closer, as in Fig.Qd. 




FIG. 3: (Color online) GP of su(2) for the largest energy 
eigenstate as a function of normalized self-collision parameter 
and the number of particles, with cross-collision rate 77 = 0. 

Our results for the phase space distribution of the max- 
imum energy state are confirmed as we analyze the be- 
havior of the GP as a function of the self-collision param- 
eter and the number of particles, as shown in Fig.(|3|, ne- 
glecting the cross-collisions. The GP initially decreases 
slowly with independent of the value of N, corre- 
sponding to the region where the distributions broadens 
along the x-axis. However, close to ^ = |, the GP 
begins to decrease quickly and, although smoothly, the 
lowering of its value is more and more steep as we increase 
N. This behavior of V su {2) suggests us a strong depen- 
dence between the derivative of the GP with respect to 
^ and the number of particles. In Fig.Q we show the 
derivative of the GP with respect to the normalized self- 
collision parameter ^ for various values of N . For an 
increasing number of particles, we see the minimum value 
of the derivative of V su {2) to move to the left side, closer 
to the classical critical value (0.5), and also the mini- 
mum becomes more pronounced. We define the value of 
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the GP for this model is an excellent measure for both 
to indicate the dynamical loss of coherence of initially 
coherent state and to indicate the QPT. When a state 
suffers fundamental changes resulting from the QPT, its 
GP must follow its behavior, because it has all the infor- 
mation about its coherence and the degree of localization 
over the phase space. 

We acknowledge J. Vidal for bringing aspects of the 
LMG model to our knowledge. Our work is supported 
by FAPESP, and CNPq. 



FIG. 4: (Color online) Derivative of the GP of the eigenstate 
of largest energy eigenvalue with respect to ^ for various 
values of the total number of particles N. The dashed ver- 
tical line represents the position at the value of the classical 
transition point — |. All the curves for finite systems 
cross exactly at this critical value. 



^ at the minimum of the derivative of GP as the criti- 
cal value of the quantum dynamical transition k%(N). It 
is already clear that the value of k%(N) is brought closer 
to the classical transition value k c = ^ for increasing 
N, but we still need to characterize how this approxi- 



mat ion happens. The values of In 



from the 



curves in Fig. Q suggest a power law between (k% — k c ) 
and N. A linear interpolation of the data points gives 

K Q _ K — il e 0,31±0,05^Y-0,657±0,009 ^ jy-l,657±0,009 jj. 

c c N 



is evident that k% 



as N — > oo. 



In conclusion we considered the GP to analyze the dy- 
namics of coherence loss of initially coherent state for the 
BECSDW model. In the MST regime the GP remains 
high without significant decoherence. For the JO regime, 
on the other hand, we have seen that the first decay of 
the GP was in a similar time scale to the MST one, be- 
ing the decoherence much more significant. In the JO 
regime no considerable re-coherence can be achieved and 
the GP correspondingly has a much lower value than in 
the previous regime at this time. Since the coherent state 
represents the closest to the classical state, the value of 
GP enabled us to estimate the quality of the semiclassical 
approximation at each time in both regimes. Moreover 
we have employed the GP as a tool for characterizing 
a QPT in the same model. We have shown for finite 
number of particles (N < oo) the bifurcation of the Q- 
function of the largest energy state as the self-collision 
parameter n becomes larger than a critical value k%{N). 
Also, we have shown the suppression of cross-collisions 
between the particles in different wells. By increasing 
the number of particles N, the GP has shown a more 
and more steeper behavior near the critical value of ^; 
moreover, its value has tended to the known classical 
value = | as N — > oo. Finally we have shown that 
n q c — k c is consistent with a power law in N. Therefore, 
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